Spectroscopic source redshifts and parameter constraints from weak lensing and CMB 
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Weak lensing is a potentially robust and model-independent cosmological probe, but its accuracy 
is dependent on knowledge of the redshift distribution of the source galaxies used. The most robust 
way to determine the redshift distribution is via spectroscopy of a subsample of the source galaxies. 
We forecast constraints from combining CMB anisotropics with cosmic shear using a spectroscopi- 
cally determined redshift distribution, varying the number of spectra N avec obtained from 64 to oo. 
The source redshift distribution is expanded in a Fourier series, and the amplitudes of each mode are 
considered as parameters to be constrained via both the spectroscopic and weak lensing data. We 
assume independent source redshifts, and consider in what circumstances this is a good approxima- 
tion (the sources are clustered and for narrow spectroscopic surveys with many objects this results 
in the redshifts being correlated). It is found that for the surveys considered and for a prior of 0.04 
on the calibration parameters, the addition of redshift information make significant improvements 
on the constraints on the cosmological parameters; however, beyond N spec ~fewxl0 3 the addition 
of further spectra will make only a very small improvement to the cosmological parameters. We 
find that a better calibration makes large N spec more useful. Using an eigenvector analysis, we find 
that the improvement continues with even higher N apec , but not in directions that dominate the 
uncertainties on the standard cosmological parameters. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 



I. INTRODUCTION 

Weak lensing (WL) is a promising tool for an era of 
precision cosmology (for reviews, see 0, 0, S an d ref- 
erences therein) . Many recent studies showed the poten- 
tial of this established technique in constraining various 
cosmological parameters 0, m B M> 1st fTot fill .112. Ha . Jl4L 
E E3. Using current data, Refs. 0Tlia 111 M EH 
showed that WL can provide constraints that are com- 
petitive with other cosmological probes. Many larger 
and more ambitious surveys are ongoing, planned or pro- 
posed. These include the Deep Lens Survey 1231 : the 
NOAO Deep Survey; the CFHT Legacy Survey [U Pan- 
STARRS; SNAP jMlHIll; and LSST HJ. A major 
challenge for the future of WL studies is to have a very 
tight control of systematic errors. These include incom- 
plete knowledge of the source redshift distribution [Hl28||. 
Recent cosmic shear studies have obtained their redshift 
distributions from joint redshift-magnitude distributions 
measured in independent surveys |2fJ; from photometric 
redshifts ("photo-z's") 29] ; or some combination. The 
former approach suffers from the difficulty that the se- 
lection function (and relative weighting) of galaxies may 
be different in the lensing and redshift surveys, while the 
latter approach suffers from possible photo- z errors that 
are difficult to constrain in the absence of spectroscopic 
confirmation, especially if only a small number of col- 
ors are measured. Several studies have marginalized over 
the source redshift distribution. Ultimately the most ro- 
bust and model-independent determination of the source 
redshift distribution would be via spectroscopy of a ran- 
domly chosen sub-sample of the source catalog. How- 
ever, spectroscopy of faint galaxies can be very time- 
consuming even with today's large- aperture telescopes. 



In this paper, we forecast Fisher-matrix constraints from 
cosmic shear using a spectroscopically determined red- 
shift distribution, varying the number of spectra N spec 
obtained. The source redshift distribution is expanded 
in a Fourier series, and the amplitudes of each mode are 
considered as parameters to be constrained via both the 
spectroscopic and WL data. We address the effect of 
galaxy clustering on our analysis and discuss the effect 
of a better calibration. 



II. MODEL PARAMETERS 

The following basic parameter set for WL is consid- 
ered: fl m h 2 , the physical matter density; J7a and w, re- 
spectively the fraction of the critical density in a dark 
energy component and its equation of state; n s (ko = 
0.05/i/Mpc) and a Sl the spectral index and running of 
the primordial scalar power spectrum at fco; £7 8 m , the 
amplitude of linear fluctuations; {<Zi---<2 , j max }, the leading 
coefficients of a Fourier expansion of the source galaxy 
redshift distribution (see Eq. 0]in the see next section); 
in the case of two bin tomography we use {qjA,QjB} 
with j = l..j m ax- We also include C, s and £ r as de- 
fined in Il 6 f to par ametrize the shear calibration bias 
@ H3, l3lll32l I33ll34| , in which the gravitational shear is 
systematically over- or under-estimated by a multiplica- 
tive factor, i.e. P K (£) = (l + ( s )P K (l), where P K (l) is the 
convergence power spectrum obtained in the absence of 
calibration errors. ( s refers to the calibration error of the 
power spectrum, which is twice the calibration error of 
the amplitude as the power spectrum is proportional to 
amplitude squared. When we consider tomography, we 
must also consider the relative calibration £ r between the 
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two redshift bins. This error affects the measured power 
spectrum P K (£) in accordance with: 



p aa = (i + / B C r )C(4 



AA i 



P, AB = (1 



Ib - Sa 



Cr)P AB (() 



p bb = 



BB i 



(1) 



where /a and /b are the fraction of the source galaxies 
in bin A and B respectively. 

In order to combine this with information from the 
CMB we include flbh 2 , the physical baryon density; t, 
the optical depth to reionization; and T/S, the tensor- 
to-scalar fluctuation ratio. We assume a spatially flat 
Universe with Q m + Q^ = 1, thereby fixing fl m and Hq as 
functions of our basic parameters, and we do not include 
massive neutrinos, or primordial isocurvature perturba- 
tions. We use as fiducial model (e.g. Ref. [35|], with w 
and T/S added): tt b h 2 = 0.0224, fl m h 2 = 0.135, f2 A = 
0.73, w = -1.0, n s = 0.93, a s = 0.0, er 8 = 0.84, r = 0.17, 
T/S = 0.2, Cs = 0.0, Cr = 0.0, and q t , q jA , q jB = 0.0. 

We will consider surveys with f s f. y — 0.01 and 0.1, 
and have number density-to-shape noise ratio n/ ("f 2 nt ) = 
4.1 x 10 9 sr , corresponding to a number density of 
n = 30 galaxies/arcmin 2 and shape+measurement noise 
(lint) = (0.3) 2 . (Note that this is the shape noise in 
the shear 7, rather than the ellipticity which is roughly 
e ss 27 if isophotal or adaptive ellipticities are used 

mmo 



III. PARAMETRIZATION OF THE REDSHIFT 
DISTRIBUTION 

The source averaged distance ratio appears as a wei ght- 
ing function in the convergence power spectrum |37l l38t 
l39| and is given by 



where n(x(z)) is the normalized source redshift distribu- 
tion. Deviations from the fiducial distribution, psl ]. 



The {qj} are thus the coefficients in the Fourier expansion 
of n(z)/n (z) in the interval < Po < 1. We have used 
the cumulative fiducial probability Po(z) instead of z as 
the independent variable because this will result in un- 
corrected constraints on the qj from spectroscopy; the 
j = coefficient vanishes because of the normalization 
constraint J n(z)dz — J no(z)dz = 1. The selection of 
the cosines is arbitrary and any complete set of orthog- 
onal functions would work just as well. The distribution 
is completely specified by the {qj}j^x, although in this 
paper we cut off the series at some j max (we show results 
for 5 and 100). For tomography, the normalized distri- 
butions and the respective Fourier expansions are given 
as above, except for bin A we replace no (z) with 

n Q (z) 



1 -5/e 2 
<{z) = 



for 



< 2z , 



for z p > 2z , 



and for bin B, we replace no(z) with 

riQ(z) = for z p < 2zq, 



'(*) = 



np(z) 
5/e 2 



for z p > 2zq. 



(0) 



(7) 



The cumulative probabilities for the two bins are 



Pfa) = fo<(z'W aud P B (z) = f*n B (z>)dz'. 



We 



use the parameters {qi..q$,q^..q^,q B ..q B } to vary the 
redshift distribution, see Fig. ^ 




n(z) = n (z) 



1 + -^=^2 qj cos(jTrP {z)) 



(3) 



(which peaks at z p = 2zq = 0.70 and has z me d w 0.94 
are parameterized via a Fourier series, 



(4) 



where Po is the cumulative redshift distribution in the 
fiducial model, 

P (z) = l-(l+- + £) e-*/*° = f n (z')dz'. (5) 



FIG. 1: Variations of the redshift distribution. The fiducial 
distribution no(z) (which peaks at z p — 2zo) and the Fourier 
expansion terms are plotted. The distributions ni±(z) rep- 
resent qi varied by ±25%. We also plot 722+ (z) to 725+ (z). 
While we varied the q'jS by ±0.05 in the analysis, we plot the 
25% variations just for the clarity of the plot. 



IV. FISHER-MATRIX ANALYSIS 



by 



The statistical error on a given parameter p a is given 



a 2 (p a ) w [(Fcmb + F WL ± F spec + n)'f a , (8) 
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where II is the prior curvature matrix, and Fcmb, Fwl 
and F spec are the Fisher matrices from CMB, WL, and 
spectroscopy, respectively. We use for F\yl the approach 
described in Rcf. with a cutoff at £ max = 3000 since 
on smaller scales the trispectrum contribution to the WL 
covariance (neglected in Eq. [BJ dominates [Zfll l4lj | . For 
CMB, we use the 4 year WMAP parameter constraints 
including TT, TE, and EE power spectra, assuming 
f sky = 0.768 (the KpO mask of Ref. H3), tempera- 
ture noise of 400, 480, and 580 /iK arcmin in Q, V, and 
W bands respectively (the rms noise was multiplied by 
\/2 for polarization), and the beam transfer functions of 
Ref. [43. The spectroscopy Fisher matrix is obtained as 
follows. If a galaxy is chosen at random from the lens- 
ing catalog, and is spectroscopically determined to have 
redshift z, then the log- likelihood for the redshift distri- 
bution parameters qj is 

\nC{qj) = lnn (z) + In 1 + ^ 1j cos 0' 7r - p o(2)) 

j 

The contribution to the Fisher matrix from this single 
galaxy is 



d\n£d\n£ 

dqj dq k 
cos(jnP (z)) cos(kirP (z)) 



%io) 



The last equality follows from the orthonormality of the 
Fourier modes, combined with the fact that Pq(z) is 
uniformly distributed between and 1. If N spec galax- 
ies have their spectra measured, and these galaxies are 
drawn independently from the source catalog, then the 
above Fisher matrix is multiplied by N spec : 



p(N spec )( cu) = -N 81. 

± spec XHji'lk! ^ ly spec u jk- 



(11) 



In the case of tomography, we have only N spec /2 spectra 
in each of the two tomography bins A and B, and so: 



F (N apec )( A M _ }_ N g 
r spec Wj ' Ik I — „ iv spec"jfc, 



(12) 



and similarly for bin B. In principle it is possible to take 
different numbers of spectra in the two bins; we have 
not attempted any optimization of this. We consider the 



cases of N s . 



0, 64, 512 and 4096 respectively. The 



results can be compared to the case where the source 
redshift distribution is known exactly by taking the limit 

Nspec * OO. 

The Fisher matrix is an asymptotic expansion, and 
while it is the standard tool of parameter forecasting, it 
sometimes leads to over-optimistic parameter constraints 
(see Ref. f° r an extreme example). We have not 
tested its validity when used with this many parameters, 
although this could be done by Monte Carlo methods as 
used by Ref. 0. 



V. SPECTROSCOPY 

In the discussion above, it has been assumed that the 
spectroscopically obtained n(z) is on average the same 
as the redshift distribution n(z) of the sources used for 
lensing, and that the redshifts of the galaxies targeted for 
spectroscopy are independent. These assumptions will be 
literally true in the idealized case that the spectroscopic 
galaxies are chosen independently from the lensing source 
catalog, and there are no failures to obtain spectro-z's. 
In practice, spectro-zs would probably be obtained by a 
multi-object spectrograph attached to a large telescope. 
The spectro-z failure rate will be minimized by using the 
longest practical integration time to maximize signal-to- 
noise for each spectrum, and the demands on telescope 
time are thus reduced if multiple objects in the same 
field of view can be targeted simultaneously. However in 
this case the galaxies are not being drawn independently 
from the source catalog, in particular large-scale cluster- 
ing can cause the redshifts of neighboring source galaxies 
to become correlated. The feasibility of obtaining many 
independent spectro-z's is directly tied to the clustering 
of the sources and the field of view of the telescope since 
these determine the maximum number of galaxies that 
can be simultaneously targeted for spectroscopy without 
the results being strongly correlated. 

The effect of source clustering can be understood 
within the context of Fisher matrix theory as follows. 
If the fiducial model is correct, and the analysis is done 
assuming that the spectro-zs are independent, the differ- 
ence Sp a = p a (est) — p"(fid) between the estimated and 
fiducial model parameters is roughly 



5p a 



<91n£ 



dpi 3 



(13) 



lid 



where C is the likelihood assuming independent spectro- 
zs (i.e. the sum of Eq. [!|] over the spectroscopic tar- 
gets), F is the Fisher matrix (also assuming indepen- 
dent spectro-zs), and the gradient is taken at the fiducial 
model. So long as each galaxy in the lensing source cat- 
alog has an equal probability of being targeted for spec- 
troscopy, it is easy to see that the spectroscopy contri- 
bution to = ^E^r(cos(j7rP (^))> = and 

hence (Sp a ) = 0, i.e. the inclusion of the spectroscopy 
likelihood function introduces no bias in the parameters 
even if the spectroscopic galaxies are not chosen inde- 
pendently. However, the clustering of the sources does 
increase the covariance matrix of the Sp a ; taking the co- 
variance of Eq. ifH^ gives 



{8p a 8p^ 



[ F -i]«7[ F - 



-11/35 



d In L d In L 



dp"< dp s 1 fid 
= [F -1 ]" 7 [F -1 ]^ 5 j [Fcmb + F WL + tt\ jS 



a,b=l 
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If we included only the a = b terms in Eq. 114(1. the 
last term in {} would simply be the spectroscopy Fisher 
matrix for independent spectro-zs, F spec ^s. The a =/= b 
terms are only non-zero due to correlations of the galaxies 
with each other, and they contribute to quantity in {} by 
an amount A g . gfe . Since the probability of two galaxies 
separated by angle 6 a b being physically associated with 
each other is oj(6 a b)/[l + u>{6 a b)], where u> is the angular 
correlation function, and the cosine is bounded in the 
range —1 to +1, it follows that the a ^ b contribution is 



WL surveys. For N s . 



512 and a field of view of 0.5 de- 



I A, 



< 



2 ^ 1 



Uj(Oab) 



a^b 



+ Uj(6 ab ) 



< 



(15) 



a^b 



as compared with F spec ^ jqk — ^N spec 6jk- Therefore if 
Eo//i u (M ^ ^spec/2, then the clustering contribution 
Ays will be small compared to the spectroscopy Fisher 
matrix F spec/y s and hence will contribute negligibly to 
the parameter uncertainties according to Eq. (|14|) . 

Let us consider the implications of this result for an 
idealized spectroscopic survey strategy that consists of 
observing J widely separated fields, and obtaining M = 
N spec /J spectra in each field. We suppose the telescope 
has a circular field of view of angular radius r, that the 
number density of galaxies targeted for selection in each 
field is n = M/irr 2 , and that the angular correlation 
function is uj(6) = (6 /6) - 7 . We then have 



a^b 




0.7 



|01 - 2 



13.25 n J J9£ 7 r- 



0.7 



where the integrals are taken over a circular disk of radius 
r. (The replacement of the galaxy-galaxy pair summa- 
tion by an integral is appropriate because the integral 
is convergent at small separation |6q — 6*2 1 • ) The fac- 
tor of J in the first line occurs because we have to re- 
peat the sum for each of the J fields. In order to have 
J2 a -ib ^{Qab) "C N spec /2, Eq. (|Po)l tells us that we need a 
number of fields given by 



./ - 2MN spec ( -i 



0.7 



(17) 



Not surprisingly, the number of fields that need to be 
observed depends on the ratio of the field-of-view ra- 
dius r to the angular clustering scale 9q. An angular 
clustering scale of 8q ~ 0.0002 degrees is observed in 
SDSS 46] for the magnitude range 21 < r < 22 (valid 
at separations of 1-30 arcmin). Even smaller 9q ap- 
plies to fainter samples, e.g. the CFDF survey ^3] finds 
w(l') w 0.01 for 18.5 < Iab < 25 galaxies, correspond- 
ing to 8 = 2.3 x 10~ 5 degrees (for slope —0.7, which 
is consistent with the u>(9) data from < 1 to several ar- 
cminutes). The fainter CFDF sample is probably more 
representative of the galaxies that will be used in future 



gree radius, we find from Eq. I|17f) the requirement J>6 
(0 O = 0.0002 deg) or J > 1.3 (0 O = 2.3 x 10~ 5 deg). In 
either case the number of fields that must be observed to 
achieve N spec — 512 ranges from a few to a few dozen, 
and for N spec — 4096 we find that the minimum number 
of fields J is a few dozen to a few hundred. The number of 
spectra to be obtained per field is M = N spec /J <~ 100. 

The problem of spectro-z failures is more difficult to as- 
sess than the source clustering. In order to accurately re- 
produce the n(z) of the lensing source catalog, the targets 
in the associated spectroscopic survey must be selected 
randomly from the lensing catalog, or at least have the 
same selection criteria. Spectro-z failures are not part of 
the lensing catalog selection criteria, and hence can bias 
the n(z) determination. The spectro-z failure rate can 
be reduced by using large- aperture telescopes and very 
long integration times, which may prove feasible if r is 
large, 8q is small, and hence the number of fields J to be 
observed is only a few. It may also be possible to reduce 
failures by imposing appropriate color cuts on the source 
sample. Nevertheless, the failure rate will never be ex- 
actly zero. The treatment of systematic errors in n(z) 
due to these failures is beyond the scope of this paper, 
but clearly deserves consideration in future work. 



VI. RESULTS 
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FIG. 2: The 68.3% confidence ellipses (assuming Fisher er- 
rors) for CMB and weak lensing. Note the improvement in the 
constraints when spectra are available, versus the case where 
we marginalize over the redshift distribution parameters. 

Our results are summarized in Tables [I] and ^ and 
Figs. [21 and |3J As expected, we find that combination of 
constraints from WL and from CMB leads to significant 
improvements in parameter estimation, notably for cs, 
f2 m /i 2 , n s , a s , w and Vl^. 

As shown in Table the increase of the number of 
expansion terms from j max — 5 to 100 has little effect on 
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TABLE I: Parameter estimation errors for WL+CMB: For WL, no-tomography, f sky — 0.01, 0.1, and ^ max = 3000. We use 
N spec = 0, 64, 512, 4096 and the limit N spec — > oo, and a priors of 0.04 on £ s and £ r (unless indicated fixed). In the last part 
of the table, we show the results with the number of terms in the series increased from 5 to a 100. 



fsky limt 1 fl b h 2 Q A <y% n„ a B T T/S w qi 92 93 ?4 95 <?((s) ^(Cr) 

CMB only-> 0.0118 0.0013 0.1981 0.086 0.059 0.036 0.018 0.187 0.872 - - - - - -__ 

AT spec= o Om 0.0096 0.0009 0.1320 0.053 0.033 0.017 0.018 0.142 0.514 4.342 19.912 34.086 66.902 58.275 0.040 0.040 

0.10 0.0092 0.0009 0.0895 0.039 0.027 0.012 0.018 0.129 0.373 1.595 6.482 13.804 21.929 18.556 0.040 0.040 

Ar spcc =64 OTOl 0.0083 0.0008 0.0248 O022 O022 CL0T2 0T017 0~T20 0~226 07162 0~244 0250 0249 O250 0~040 0.040 

0.10 0.0075 0.0007 0.0141 0.018 0.019 0.010 0.017 0.114 0.186 0.121 0.218 0.249 0.248 0.250 0.040 0.039 

7V spcc =512 Pi 0.0076 0.0008 0.0219 0.021 0.021 0.011 0.016 0.115 0.220 0.081 O088 O088 O088 O088 0.040 0.039 

0.10 0.0058 0.0007 0.0116 0.014 0.016 0.010 0.015 0.098 0.140 0.074 0.087 0.088 0.088 0.088 0.039 0.038 

Af spcc= 4096 0.01 0.0075 0.0007 0.0211 0.021 0.021 0.011 0.016 0.114 0.219 0.031 0.031 0.031 0.031 0.031 0.040 0.039 

0.10 0.0049 0.0006 0.0111 0.013 0.014 0.009 0.015 0.090 0.118 0.030 0.031 0.031 0.031 0.031 0.039 0.038 

N sp ec 00 0.01 0.0074 0.0007 0.0210 0.021 0.020 0.011 0.016 0.113 0.219 - 0.040 0.039 

0.10 0.0047 0.0006 0.0110 0.012 0.013 0.009 0.015 0.088 0.113 - - - 0.039 0.038 

N 3pcc =0 OM 0.0096 0.0009 0.1320 00511 0033 0.017 0.018 0142 0.514 4.341 19.912 34.086 66.902 58.275 - - 



Ca,r 


fixed 


0.10 


0.0092 


0.0009 


0.0895 


0.039 


0.027 


0.012 


0.018 


0.129 


0.373 


1.594 


6.482 


13.804 


21.928 


18.555 


N 3p 

<s,r 


c=64 
fixed 


0.01 
0.10 


0.0083 
0.0075 


0.0008 
0.0007 


0.0243 
0.0137 


0.022 
0.018 


0.022 
0.019 


0.012 
0.010 


0.017 
0.017 


0.120 
0.114 


0.226 
0.185 


0.159 
0.114 


0.244 
0.217 


0.250 
0.249 


0.249 
0.248 


0.250 
0.250 


N spc 

Ca,r 


c=512 
fixed 


0.01 
0.10 


0.0076 
0.0056 


0.0008 
0.0007 


0.0211 
0.0110 


0.021 
0.014 


0.021 
0.015 


0.011 
0.009 


0.016 
0.015 


0.115 
0.096 


0.220 
0.137 


0.081 
0.073 


0.088 
0.087 


0.088 
0.088 


0.088 
0.088 


0.088 
0.088 


AT spec= 4096 
£ Sjr fixed 


0.01 
0.10 


0.0074 
0.0045 


0.0007 
0.0006 


0.0202 
0.0105 


0.021 
0.012 


0.020 
0.013 


0.011 
0.009 


0.016 
0.014 


0.113 
0.085 


0.219 
0.110 


0.031 
0.030 


0.031 
0.031 


0.031 
0.031 


0.031 
0.031 


0.031 
0.031 



N spec -» 00 0.01 0.0073 0.0007 0.0201 0.021 0.020 0.011 0.016 0.113 0.218 - - - - - 
C 3 , r fixed 0.10 0.0042 0.0006 0.0104 0.012 0.012 0.009 0.014 0.083 0.104 - - - - - - 

f,ky limt' Slfr/t' 2 Q A era "s a s t T/S w qi 92 to 999 9100 o-(Cs) o-(Cr) 



N 3pi!i: 


=64 


0.01 
0.10 


0.0083 
0.0075 


0.0008 
0.0007 


0.0249 
0.0144 


0.022 
0.018 


0.022 
0.019 


0.012 
0.010 


0.017 
0.017 


0.120 
0.114 


0.227 
0.186 


0.162 
0.122 


0.244 
0.218 


0.250 
0.250 


0.250 
0.250 


0.040 
0.040 


0.040 
0.039 


N 3 pec = 


=512 


0.01 
0.10 


0.0076 
0.0058 


0.0008 
0.0007 


0.0219 
0.0116 


0.021 
0.014 


0.021 
0.016 


0.011 
0.010 


0.016 
0.015 


0.115 
0.098 


0.220 
0.140 


0.081 
0.074 


0.088 
0.087 


0.088 
0.088 


0.088 
0.088 


0.040 
0.039 


0.039 
0.038 


N S pec= 


4096 


0.01 
0.10 


0.0075 
0.0049 


0.0007 
0.0006 


0.0211 
0.0111 


0.021 
0.013 


0.021 
0.014 


0.011 
0.009 


0.016 
0.015 


0.114 
0.090 


0.219 
0.118 


0.031 
0.030 


0.031 
0.031 


0.031 
0.031 


0.031 
0.031 


0.040 
0.039 


0.039 
0.038 


/V 

Jv spec 


00 


0.01 
0.10 


0.0074 
0.0047 


0.0007 
0.0006 


0.0210 
0.0110 


0.021 
0.012 


0.020 
0.013 


0.011 
0.009 


0.016 
0.015 


0.113 
0.088 


0.219 
0.113 










0.040 
0.039 


0.039 
0.038 





0.9 




FIG. 3: Same as Fig. [5] but including tomography. 

this result. This suggests that the parameter estimates 
have converged and, since Eq. (@J can model an arbitrary 
function for sufficiently large j m ax, it suggests that the 
form (Eq.0} for the redshift distribution is not artificially 
constraining the cosmological parameters. 



VII. DISCUSSION 

In agreement with previous results, we find that the 
addition of redshift information is helpful for several of 



TABLE II: Parameter estimation errors for WL+CMB: WL 
fsky is varied from 0.01 to 0.1, and ^ max = 3000. We use 
N 3pec = 0, 64, 512, 4096 and the limit N spe c — * 00, and priors 
of 0.04 on £ s and ( r (unless indicated fixed). Tomgraphy case. 





fsky 


O a 


CT8 


w 


CMB only-> 




0.1981 


0.086 


0.872 


■'V spec — 


0.01 


0.0785 


0.035 


0.339 




0.10 


0.0340 


0.024 


0.240 


N spcc = 64 


0.01 


0.0169 


0.018 


0.185 




0.10 


0.0102 


0.013 


0.118 


N Bp ec = 512 


0.01 


0.0160 


0.017 


0.173 




0.10 


0.0084 


0.010 


0.088 


Ns P ec = 4096 


0.01 


0.0154 


0.016 


0.162 




0.10 


0.0078 


0.009 


0.074 


N Bp ec — * OO 


0.01 


0.0148 


0.015 


0.154 




0.10 


0.0071 


0.008 


0.066 


Nspec = 


0.01 


0.0684 


0.032 


0.311 


Cs,r fixed 


0.10 


0.0288 


0.023 


0.230 


Nspec = 64 


0.01 


0.0147 


0.017 


0.168 


Q S) r fixed 


0.10 


0.0091 


0.012 


0.108 


Nspec - 512 


0.01 


0.0137 


0.016 


0.153 


Q Sj7 . fixed 


0.10 


0.0069 


0.009 


0.075 


Nspec = 4096 


0.01 


0.0130 


0.014 


0.140 


Q Sj7 . fixed 


0.10 


0.0062 


0.008 


0.063 


Nspec — > 00 


0.01 


0.0118 


0.012 


0.122 


Qs !r fixed 


0.10 


0.0043 


0.005 


0.043 



the cosmological parameters such as Oa, n s , and w. 
The optical depth t is less improved by lensing because 
most of the statistical power on r is coming from the 
CMB polarization reionization peak, which is not degen- 
erate with any lensing-related quantities. The baryon 
density is well-constrained by the CMB, however by im- 
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proving our constraints on several cosmological param- 
eters, lensing information breaks the relatively weak re- 
maining degeneracies in the CMB and reduces the error 
from a(n b h 2 ) = 0.0013 to 0.0009. Even lensing with 
N sp ec = is sufficient to break this degeneracy, so inclu- 
sion of redshift information adds little for Q,bh 2 . 

The addition of redshift information leads to further 
significant improvements in the parameter estimation as 
expected. The N spec = constraints are substantially 



worse than N spec = 64, because in the N s . 



case 



even wildly oscillating redshift distributions n(z) are for- 
mally allowed, as is evidenced by the large a(qj) ^> 1 
in Table U Although the gradual addition of spectro- 
scopic redshifts does make additional improvements on 
the constraints on the cosmological parameters, we find 
that there is a certain number of N spec (several thou- 
sands for the surveys considered here) beyond which the 
addition of further spectra will make only a very small 
improvement to the cosmological parameters. 

In order to try to explain this, let us recall that error 
bars are correlated, and hence there can still be com- 
binations of parameters that are degraded by incom- 
plete knowledge of the source redshift distribution. In- 
formation on this degradation is contained within the 
9x9 covariance matrix C of the cosmological parame- 
ters {Q m h 2 , ttbh 2 , Qaj 08, n s> a s' T > T/S, w }- We examine 
the degradation factor R of the combination of parame- 
ters x = k a p a given by 



R(x) 



a 2 (x;N spec ) 



a 2 {x, N spec = oo) [C(oo)] a ^k a k 







(18) 



The maximum value of R can be found by setting 
dR/dk a — 0, which leads to the eigenvalue equation 

[w'' 1 - " [ c "<»>MC(AW)rt. 

(M) 

Thus the k a that maximizes the degradation R(x) is an 
eigenvector of [C _1 (oo)C(iV S p ec )], and the degradation 
factor R(x) is the eigenvalue. (The maximum degrada- 
tion factor corresponds to the maximum eigenvalue of 
[C~ 1 (oo))C(-/V spec )]; the other eigenvectors are minima 
or saddle points of R.) For N spec = 512, f s f. y = 0.1, 
and no-tomography, we find a maximum eigenvalue of 
X ma x = 2.2 with the corresponding combination of pa- 
rameters 

x = -4.37fi m /i 2 + 8.89ft A - 5.39cr 8 - 2.57n s + 4.32a s 
+0.53w + 3.lin b h 2 - 4.59r - 7.18T/S. (20) 

Having only a finite number of spectra is degrading the la 
WL+CMB constraint on x by a factor of \fX ma x » 1.5, 
but from Table [Q wo can sec that the constraints on the 
standard set of cosmological parameters is degraded by 
< 10%. This is because x is a direction in parameter 
space that is very well-constrained by WL+CMB: with 
CMB only, we have <j(x) = 1.02, whereas with WL and 
N spec = 512 added to the CMB we have a(x) = 0.098 and 



correlation coefficient p(x,qi) — —0.64. Consequently, 
although x is degraded by imperfect knowledge of the 
source redshift distribution, and is degenerate with the 
redshift distribution parameters, it is not the direction 
that dominates the uncertainties on the individual cos- 
mological parameters. The second-largest eigenvalue of 
[C _1 (oo))C(iV ;i p ec )] is 1.2, indicating that the directions 
other than x are not significantly degraded. When we 
fix the calibration parameters, the eigenvector analysis 
shows an increase of almost an order of magnitude on 
the degradation with y/ X max w 16.6, but again, not in 
the direction that dominates the uncertainties on the con- 
ventional cosmological parameters. As seen from Table 
HI and Table |n] a better calibration will make a larger 
numbers of spectra slightly more useful but a significant 
improvement require a rather idealized perfect knowledge 
of the calibration. 

The most significant improvement in the cosmological 
parameters with large numbers of spectra comes from the 
cases with tomography, since this helps break the degen- 
eracies that dominate the parameter uncertainties and 
allows the directions well constrained by WL to become 
more important. One can see from Table [H] that the 
uncertainties on ag, and w are degraded by ~ 50% 



for N. 



spec 



4096 versus N. 



spec 



sky 



oo with the larger area 
0.1. In all the other cases the degradation with 
N spec = 4096 is much less. 

It is possible that different results would be obtained if 
WL were combined with additional cosmological probes 
such as supernovae, which could break remaining degen- 
eracies in the data and therefore make the degradation 
of x more important. A general way to investigate this 
possibility is to examine X max for the WL-only matrices, 
since combining WL with cosmological probes that do 
not include the {qj} must result in a degradation factor 
X ma x smaller than that for WL alone. When we apply 
the eigenvalue analysis to the weak lensing only matrices, 
we obtain the degradation factors of \/X ma x ~ 19.98 for 
the fixed calibration case and N spec = 512. The lat- 
ter number indicates that there are directions in cos- 
mological parameter space which are dramatically im- 
proved by exact knowledge of the redshift distribution 
rather than only 512 spectra. We can reduce the degra- 
dation factor to X max = 2 by increasing N spec to 200000, 
and it is thus only for N spec > 200000 that every di- 
rection in parameter space is limited by lensing statis- 
tics rather than redshift distribution uncertainties. How- 
ever, this improvement is rapidly lost due to the approx- 
imate calibration-rcdshift degeneracy if the calibration is 
uncertain: if the prior on the calibration parameters is 
widened to cr(Cs) = cr(Cr) = 0.02 we achieve X ma x = 2 
at N spec = 4000, and with a(( s ) = a(( r ) = 0.04 we 
achieve X ma x — 2 at N spec = 1100. Thus having either 
good calibration or a well-determined redshift distribu- 
tion individually may not be very useful, but having both 
combined can significantly improve some constraints. 

We conclude that, though significant improvement is 
obtained from the addition of redshifts information, there 



7 



is a certain number of N spec - of order 10 3 for the surveys 
considered here - beyond which the addition of further 
spectra will make only a very small improvement to the 
cosmological parameters. We do find that there are direc- 
tions in parameter space that continue to improve and do 
not saturate until N spec is very large, especially if the cali- 
bration is very well determined; however these directions 
correspond to very large eigenvalues of the WL Fisher 
matrix and do not dominate the parameter uncertainties 
in the WL+CMB combinations we have discussed. The 
results presented here indicate that if N spec ~fewxl0 3 
spectra of representative sources can be obtained, cos- 
mological constraints can potentially be improved rela- 
tive to the CMB-only case without model-dependent as- 
sumptions about the source redshift distribution. These 
results are robust against fluctuations in the spectro-z 



distribution due to galaxy clustering so long as enough 
spectroscopic fields are observed, as detailed in Sec. [3 
while future work is required to address the spectro-z 
failures. 
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